Search results for "Computer Science::Numerical Analysis"
showing 10 items of 32 documents
An algebraic multigrid based shifted-Laplacian preconditioner for the Helmholtz equation
2007
A preconditioner defined by an algebraic multigrid cycle for a damped Helmholtz operator is proposed for the Helmholtz equation. This approach is well suited for acoustic scattering problems in complicated computational domains and with varying material properties. The spectral properties of the preconditioned systems and the convergence of the GMRES method are studied with linear, quadratic, and cubic finite element discretizations. Numerical experiments are performed with two-dimensional problems describing acoustic scattering in a cross-section of a car cabin and in a layered medium. Asymptotically the number of iterations grows linearly with respect to the frequency while for lower freq…
A damping preconditioner for time-harmonic wave equations in fluid and elastic material
2009
A physical damping is considered as a preconditioning technique for acoustic and elastic wave scattering. The earlier preconditioners for the Helmholtz equation are generalized for elastic materials and three-dimensional domains. An algebraic multigrid method is used in approximating the inverse of damped operators. Several numerical experiments demonstrate the behavior of the method in complicated two-dimensional and three-dimensional domains. peerReviewed
Wulff shape characterizations in overdetermined anisotropic elliptic problems
2017
We study some overdetermined problems for possibly anisotropic degenerate elliptic PDEs, including the well-known Serrin's overdetermined problem, and we prove the corresponding Wulff shape characterizations by using some integral identities and just one pointwise inequality. Our techniques provide a somehow unified approach to this variety of problems.
On a global superconvergence of the gradient of linear triangular elements
1987
Abstract We study a simple superconvergent scheme which recovers the gradient when solving a second-order elliptic problem in the plane by the usual linear elements. The recovered gradient globally approximates the true gradient even by one order of accuracy higher in the L 2 -norm than the piecewise constant gradient of the Ritz—Galerkin solution. A superconvergent approximation to the boundary flux is presented as well.
Functional Type Error Control for Stabilised Space-Time IgA Approximations to Parabolic Problems
2018
The paper is concerned with reliable space-time IgA schemes for parabolic initial-boundary value problems. We deduce a posteriori error estimates and investigate their applicability to space-time IgA approximations. Since the derivation is based on purely functional arguments, the estimates do not contain mesh dependent constants and are valid for any approximation from the admissible (energy) class. In particular, they imply estimates for discrete norms associated with stabilised space-time IgA approximations. Finally, we illustrate the reliability and efficiency of presented error estimates for the approximate solutions recovered with IgA techniques on a model example.
(Approximate) Low-Mode Averaging with a new Multigrid Eigensolver
2015
We present a multigrid based eigensolver for computing low-modes of the Hermitian Wilson Dirac operator. For the non-Hermitian case multigrid methods have already replaced conventional Krylov subspace solvers in many lattice QCD computations. Since the $\gamma_5$-preserving aggregation based interpolation used in our multigrid method is valid for both, the Hermitian and the non-Hermitian case, inversions of very ill-conditioned shifted systems with the Hermitian operator become feasible. This enables the use of multigrid within shift-and-invert type eigensolvers. We show numerical results from our MPI-C implementation of a Rayleigh quotient iteration with multigrid. For state-of-the-art lat…
On the proper homotopy invariance of the Tucker property
2006
A non-compact polyhedron P is Tucker if, for any compact subset K ⊂ P, the fundamental group π1(P − K) is finitely generated. The main result of this note is that a manifold which is proper homotopy equivalent to a Tucker polyhedron is Tucker. We use Poenaru’s theory of the equivalence relations forced by the singularities of a non-degenerate simplicial map.
Generalized Alomari functionals
2015
We consider a generalized form of certain integral inequalities given by Guessab, Schmeisser and Alomari. The trapezoidal, mid point, Simpson, Newton-Simpson rules are obtained as special cases. Also, inequalities for the generalized Alomari functional in terms of the $n$-th order modulus, $n=\overline{1,4}$, are given and applied to some known quadrature rules.
NUMERICAL ALGORITHMS
2013
For many systems of differential equations modeling problems in science and engineering, there are natural splittings of the right hand side into two parts, one non-stiff or mildly stiff, and the other one stiff. For such systems implicit-explicit (IMEX) integration combines an explicit scheme for the non-stiff part with an implicit scheme for the stiff part. In a recent series of papers two of the authors (Sandu and Zhang) have developed IMEX GLMs, a family of implicit-explicit schemes based on general linear methods. It has been shown that, due to their high stage order, IMEX GLMs require no additional coupling order conditions, and are not marred by order reduction. This work develops a …
Stability of radial symmetry for a Monge-Ampère overdetermined problem
2008
Recently the symmetry of solutions to overdetermined problems has been established for the class of Hessian operators, including the Monge-Ampère operator. In this paper we prove that the radial symmetry of the domain and of the solution to an overdetermined Dirichlet problem for the Monge-Ampère equation is stable under suitable perturbations of the data. © 2008 Fondazione Annali di Matematica Pura ed Applicata and Springer-Verlag.